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Abstract. We consider evolution in the unit disk in which the sample paths are rep- 
resented by the trajectories of points evolving randomly under the generalized Loewner 
equation. The driving mechanism differs from the SLE evolution, but nevertheless solu- 
tions possess similar invariance properties. 



1. Introduction 

The last decade has been marked by a burst of interest in Schramm- Lowner Evolution 
(SLE). SLE has led to an elegant description of several 2D conformally invariant statistical 
physical systems at criticality by means of a solution to the initial value problem for a 
special differential equation with a random driving term given by ID Brownian motion. 
The origin of SLE can be traced to the seminal Schramm's paper [22], where he revisited 
the notion of scaling limit and conformal invariance for the loop erased random walk and 
the uniform spanning tree. This approach was developed in many further works and many 
lattice physical processes were proved to converge to SLE with some specific diffusion 
factor, e.g., percolation, Ising model [HJ ESI ES]- On the other hand, SLE became an 
approach to conformal field theory which emphasizes CFT's roots in statistical physics 

We consider another setup in which the sample paths are represented by the trajectories 
of points (e.g., the origin) in the unit disk D evolving randomly under the generalized 
Lowner equation. The driving mechanism differs from SLE. In the SLE case the Denjoy- 
Wolff attracting point (oo in the chordal case or a boundary point of the unit disk in the 
radial case) is fixed. In our case, the attracting point is the driving mechanism and the 
Denjoy- Wolff point is different from it. In analytic terms, the generalized Lowner evolution 
{4>t}t>o solves the initial value problem 




(L) j, ,rM--W2)-rW)(r(^)-l)p(^) ) t), z( _ B ^ t>Q) 
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where r : [0, +00) — » O is measurable and the driving function p(z, t) is a so-called Herglotz 
function (see Section [2] for the definition). 
In the radial SLE case r = and 



In our case r(t) = e lkBt , and the Denjoy- Wolff point is different from r, and there is no 
common fixed point of the evolution. 

Even though the driving mechanism in our case differs from that of SLE, the generated 
families of conformal maps still possess the important time-homogeneous Markov property. 

In the deterministic case, Section 3, we thoroughly study a particular case when r(t) = 
e lkt , k £ M. We give a complete description of the case of the evolution of Mobius au- 
tomorphisms of D. In particular, we obtain the intervals of the constant k in which the 
evolution is elliptic, hyperbolic and parabolic. In the general case of endomorphisms of D 
we prove a sufficient condition of ellipticity of the dynamics for sufficient large k. We also 
find values of the parameter k for which the trajectories of a point are closed. 

The proposed model also describes deterministic and stochastic evolution of shapes in 
the complex plane. Following [23], we understand by a shape a simple closed smooth curve. 
The study of 2D shapes is one of the central problems in the field of computer vision. A 
program of such study and its importance was summarized by Mumford at ICM 2002 in 
Beijing [17]. However, we plan to focus on shape evolution in forthcoming papers. 

The central part of the paper, Section 4, is dedicated to a stochastic version of ([I]) in 
which the driving term is r(t) = e lkBt , where B t is the standard 1-dimensional Brownian 
motion. The equation ([T]) is reduced to an Ito 2D diffusion SDE, which we solve explicitly 
in some cases. We calculate the infinitesimal generator which is given by the Virasoro 
generators in a certain representation. We give estimates of growth of calculate 
momenta, and describe the boundary diffusion on T. 

So far, the only well-studied examples of the generalized Lowner equation were the 
classical cases corresponding to r(£) = (radial Lowner equation) and r(t) = 1 (chordal 
Lowner equation). This paper is one of the first studies where examples of generalized 
Lowner equation with non-constant attracting point r(t) are considered. 

Acknowledgments. The authors would like to thank Filippo Bracci and David Shoikhet 
for inspiring discussions. The authors are especially grateful to Alan Sola for reading the 
original version of the paper and for his important critical remarks. 



The pioneering idea of Lowner [TB] in 1923 contained two main ingredients: subordina- 
tion chains and semigroups of conformal maps. This far-reaching program was created in 
the hopes to solve the Bieberbach conjecture [I] and the final proof of this conjecture by 
de Branges [9] in 1984 was based on Lowner's parametric method. The modern form of this 
method is due to Kufarev [13] and Pommerenke [191 120] . We concentrate our attention on 



p(z,t) 




where B t is Brownian motion. 



2. Generalized Lowner evolution 
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a generalization of the Lowner-Kufarev approach by F. Bracci, M. D. Contreras, S. Diaz- 
Madrigal, and P. Gumenyuk [5J [S]. One of the main results of the generalized Lowner 
theory is an essentially one-to-one correspondence between evolution families, Herglotz 
vector fields and (generalized) Lowner chains. Below we briefly describe the basic notions 
of the theory. 

2.1. Preliminaries: Semigroups and fixed points of holomorphic functions. By 

the Schwarz-Pick lemma, every holomorphic self-map <p of the unit disk D may have at most 
one fixed point r in D. If such a point r exists, then the point is called the Denjoy- Wolff 
point of <p, and in this case the function (p is said to be an elliptic self-map of D. 

Otherwise, there exists a point r on the unit circle T, such that the angular limit 
Z lim z ^ T <p(z) = t (this statement is known as the Denjoy- Wolff theorem). The point 
t is again called the Denjoy- Wolff point of <p. The limit Z lim z _ s>T <p'(z) = ay always exists, 
its value a v belongs to (0, 1], and the map ip in this case is said to be either hyperbolic (if 
a 9 < 1) or parabolic (if a v = 1) (for details and proofs see, e. g., [IT 124"]). 

A family {ipt}t>o of holomorphic self-maps of the the unit disk D is called a one-parameter 
continuous semigroup if 

(1) tp = id D , 

(2) i) t+s = %p t oip sl for s,t>0, 

(3) lining i])t{z) = 4>s( z ), for all s > and z£D, 

(4) limf_;.o+ il>t(z) = z locally uniformly in D. 

For every semigroup {ijjt{z)}t>o there exists a unique holomorphic function / : © — > C, 
such that the semigroup {ipt(z)} t >o is the unique solution to the initial- value problem 



The function / is called the infinitesimal generator of the semigroup {ipt}t>o- 

A holomorphic function / is an infinitesimal generator of a semigroup if and only if there 
exists r G © and a holomorphic function p : D — > C with Rep > 0, such that 



This representation is unique (if / ^ 0), and is known as the Berkson- Porta representation 
of /. The point r turns out to be the Denjoy- Wolff point of all functions in {ipt( z )}t>o- We 
will use the term Herglotz function for the function p(z). 

We say that a semigroup ipt(z) is elliptic, hyperbolic or parabolic if it consists of elliptic, 
hyperbolic or parabolic self-maps of D, respectively. 

2.2. Evolution families and Herglotz vector fields. We start with the notion of evo- 
lution family. 

Definition 1. An evolution family of order d G [l,+oo] is a two-parameter family 
{0s,t}o<s<*<+oo of holomorphic self-maps of the unit disk, such that the following three 
conditions are satisfied. 



(2) 



ip (z) = z. 



f(z) = (z - t)(tz - l)p(z) 



• SjS = /'(/- : 
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• </>s,t = 0it,i ° 0s,« f or 0,11 < S < U < t < +00; 

• for any z G D and T > there is a non-negative function k Z) T G L d ([0, T], R), swc/t 

|<M^) - <M*)I < / k z>T {£)d£, z G D 

J u 

for all < s < u <t <T. 

An infinitesimal description of an evolution family is given in terms of a Herglotz vector 
field. 

Definition 2. A (generalized) Herglotz vector field of order d is a function G : D x 
[0, +00) — > C satisfying the following conditions: 

• the function [0, +00) 3 t 1— > G(z, t) is measurable for every z G D; 

• the function z 1— > G(z, t) is holomorphic in the unit disk for t G [0, +00); 

• for any compact set K C D and for every T > there exists a non-negative function 
k K ,T G L d ([0,T],R), such that 

\G(z,t)\ < k K , T {t) 

for all z G K, and for almost every t G [0, T]; 

• for almost every t G [0, +00) the vector field G(-,t) is semicomplete. 

By semicompleteness we mean that the solution to the problem 

dx{l) 



dl 

x(s) 



G(x(l),t), 



is defined for all times / G [s, +00), for any fixed s > 0, fixed t > and fixed z6D. 

An important result of general Lowner theory is the fact that the evolution families can 
be put into a one-to-one correspondence with the Herglotz vector fields by means of the 
so-called generalized Lowner ODE. This can be formulated as the following theorem. 

Theorem 1 ([5, Theorem 1.1]). For any evolution family {<p s ,t} of order d > 1 in the unit 
disk there exists an essentially unique Herglotz vector field G(z, t) of order d, such that for 
all z G © and for almost all t G [0, +00) 

(3) — j t — =G{(f> att (z),t). 

Conversely, for any Herglotz vector field G(z, t) of order d > 1 in the unit disk there exists 
a unique evolution family of order d, such that the equation above is satisfied. 

Essential uniqueness in the theorem above means that for any other Herglotz vector 
field H(z, t) satisfying d3J), the equality H(z, t) = G(z, t) holds for all z G D and almost all 
t G [0,+oo). 

Herglotz vector fields admit a convenient Berkson-Porta-like representation using so- 
called Herglotz functions. 
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Definition 3. A Herglotz function of order d G [1, +oo) is a function p : D x [0, +oo) — > C, 
such that 

• the function t \-t p(z, t) belongs to Lf oc ([0, +oo), C) for all z G D; 

• the function z i-> p(z, t) is holomorphic in D for every fixed t G [0, +oo); 

• Rep(z, £) > /or all z G © and for all t G [0, +oo) . 

Now, the representation of Herglotz vector fields is given in the following theorem. 

Theorem 2 ([5J Theorem 1.2]). Given a Herglotz vector field of order d > 1 in the unit 
disk, there exists an essentially unique (i. e., defined uniquely for almost all t for which 
G(-, t) ^ 0) measurable function r : [0, +oo) — > D and a Herglotz function p(z, t) of order 
d, such that for all z G D and almost all t G [0, +oo) 

(4) G(z,t) =(z- r{t))W)z - l)p(z,t). 

Conversely, given a measurable function r : [0, +oo) — > D and a Herglotz function p(z, t) 
of order d > 1, the vector field defined by the formula above is a Herglotz vector field of 
order d. 

According to TheoremfTJ to every evolution family {<p s ,t} one can associate an essentially 
unique Herglotz vector field G(z,t). The pair of functions (p, r) representing the vector 
field G(z,t) is called the Berkson-Porta data of the evolution family {4> s ,t\- 

Let {4> s ,t} be an evolution family with Berkson-Porta data (p,r). In the simplest case 
when neither p, nor r change in time (i. e., the corresponding Herglotz vector field G(z,t) 
is time- independent), r turns out to be precisely the Denjoy- Wolff point of every self-map 
in the family {<p s ,t}- Moreover, for any < s < +oo, we have that (f) s ,t(z) T uniformly 
on compacts subsets of D, as t — > +oo. For this reason, we call r the attracting point of 
the evolution family {<p s ,t}- 

In the case when the Herglotz field G(z,t) is time- dependent, the meaning of r is ex- 
plained in the following theorem. 

Theorem 3 ([Q Theorem 6.7]). Let {<p s ,t} be an evolution family of order d > 1 in the 
unit disk, and let G(z, t) = (z — T{t)){r{t)z — l)p(z, t) be the corresponding Herglotz vector 
field. Then for almost every s G [0, +oo) ; such that G(z,s) ^ 0, there exists a decreasing 
sequence {t n (s)} converging to s, such that $ s ,t n (s) 7^ id® and 

t(s) = lim r(s, n), 

n— >oo 

where r(s,n) denotes the Denjoy- Wolff point of 4> s ,t n (s)- 

Note, however, that it is not true in general that r(t) is the Denjoy- Wolff point of the 
map <fis,t(z) (see, for instance, Example 1 in Section [3]). 

3. Deterministic case 

We start by considering a simpler deterministic setting r(t) = e lkt . In the next Section H] 
we switch to the stochastic setting and let r(t) move as a Brownian particle on the unit 
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circle T, i.e., r(i) = e lkBt , where k > 0, and B t is a 1-dimensional standard Brownian 
motion. 

Due to the definition of an evolution family (Definition [TJ, = <f> t o and one may 
restrict attention to the one-parameter family of maps {(pt}t>o '■= {<Po,t}t>o, rather than 
consider the two-parameter evolution family. So we are going to analyze the following 
initial value problem 

U Mz) = m^r p{Mz)A ^ z£Dj (>0 

[0o(^) = 2, 

If the Herglotz function p(z,t) can be written in the form p(z,t) = p(z/r(t)), where 
p : D — 7- C is a holomorphic function with non-negative real part, then the problem ([5]) is 
simplified significantly. Indeed, the change of variables 

leads to the following initial value problem for a separable differential equation 



(6) 



so that the solution to (jSJ) can be written in quadratures. 

Observe that the vector field f(z) = (z — l) 2 p(z) — ikz in the right-hand side of ([6]) is an 
infinitesimal generator, and (EJ) is, in fact, an autonomous Lowner equation. It generates 
a semigroup {i/Jt}t>o of self- maps of the unit disk. We denote the Berkson-Porta data 
associated with this semigroup by (p (z),T ). 

If tq lies inside D, then /(to) = 0. Otherwise, if r is on the circle T, then we have 
Zlim 2 _;. T0 f(z) = 0. Hence, the corresponding Herglotz function in the Berkson-Porta de- 
composition of / is given by 

. . (z— l) 2 p(z) — ikz 
(z - t )(t z - 1) 

Example 1. The simplest example is given by the Herglotz function p{w) = which 

Tno T munor TlTiTV. nDPnmDC Iitioqt 

T(t)-Z- 



corresponds to p(z,t) = ffl. ■ The Lowner ODE becomes linear 



-j t <t>t{ z ) = T -M z )i 



and the solution to the problem (jSJ) is given by 



Uz) = e- t z+ e 



l + ik 

Let DW^(t) denote the Denjoy- Wolff point of the map 4>t- Solving the equation 4>t{z) 
we easily get that in this example 

_ ]_ _|_ gt+ikt 



DW^(t) 



(l + ik) (e* - L 
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This may be regarded as a simple example showing that in general r(t) ^ DW(t). 
The corresponding autonomous Lowner equation is 

and its solution is 

i{, t ( z ) = e ~ t{ - 1+ik h + * . 

v ' l + ik 

Let us analyze the dynamical properties of this evolution. The Denjoy- Wolff point of 
the evolution ipt( z ) is t = jh^- The point r is hyperbolic if k = 0, and elliptic otherwise. 
The Herglotz function in the Berkson-Porta decomposition of / is 

Po{z) = 



z-(l-ik)' 

Example 2. In the previous example, the solution is given by a family of Mobius transfor- 
mations (in particular, circles are always transformed into circles). If one is interested in 
the evolution of shapes, it could be more relevant to look at examples where p(z, t) is a 
transcendental analytic function with non-negative real part. Figure [1] shows the evolution 
of the domains </>t(B) corresponding to the function p(z,t) = exp ?p. 




(a) t = \ (b) t = I (c) t = f 

Figure 1. Evolution of t (D) for p(z,t) = exp^/^r), k = 1 (the point in 
T represents t(£)). 



A Java applet demonstrating examples with different Herglotz functions and different 
values of k is available at |http : / /folk . uib . no/ giv023/loew nershapes . 

3.1. Families of automorphisms. Let us describe the class of Herglotz functions gener- 
ating families of automorphisms of the unit disk, using several basic results of generation 
theory of semigroups of holomorphic self-maps of the unit disk. A summary of these results 
can be found, e. g., in [TTJ Chapter 2}. 
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Proposition 1. Let the semigroup {ipt{ z )}t>o defined by ([6]) for some k = k G R be a 
family of automorphisms of D. Then {tp t { z )}t>o is a family of automorphisms © for all 
keR. 

Proof. Note that since {ipt{z)}t>o consists of automorphisms D for k = ko, the vector field 
fi(z) = (1 — z) 2 p{z) — ik z is complete (i. e., the solution of the corresponding differential 
equation exists and is unique for all — oo < t < +oo). Also note that for a k G R the 
vector field f2(z) = i{^o — k)z is complete. Now, the vector field fz{z) = (1 — z) 2 p(z) —ikz 
is complete as a sum of two complete vector fields: fs = f\ + fi, and thus, the generated 
semigroup ?pt(z) is a family of automorphisms. □ 

Proposition 2. The family {ip t {z)}t>o consists of automorphisms if and only if the corre- 
sponding Herglotz function is of the form p(z) = A j^ 2 - + Bi, where A > 0, B G R. 

Proof. As the previous proposition suggests, it is sufficient to prove the proposition for 
k = 0. Let f(z) be the generator of {ipt(z)} t >o for k — 0, i. e., f(z) = (1 — z) 2 p(z). 
Since f(z) is a complete vector field, it must be a polynomial of the form 

f[z) = -az 2 + ibz + a, a G C, b G R, 

(see, e. g., [HI Theorem 2.6]). Since 1 is a fixed point of the semigroup, it follows from the 
relation f(l) = 0, that b = — 2Ima. Hence, f(z) is a polynomial of the form 

f(z) = {-A + Bi) z 2 - 2Biz + (A + Bi), 

which is precisely (1— z) 2 p{z) with p{z) = Aj^+Bi. The function p(z) has nonnegative real 
part if and only if A > 0. (The case A < corresponds to the Berkson- Porta decomposition 

of f(z) with the attracting point — ^zff and the Herglotz function p(z) = ^{A^iJB)+A+iB ) • 



Since p(0) = A + we will also write the Herglotz functions generating families of 
automorphisms as p(z) = y±| Rep(0) + i lmp(0). 

The differential equation for the semigroup generated by p(z) = Aj^ + Bi is written 

as 

(7) <Hh = , 

U (-A + iB) ^1 + (-2iB -ik)^j t + A + iB 

Denote by D the discriminant of the polynomial in the denominator of the left-hand 
side: D = 4A 2 — ABk - k 2 . 

D>0, if 2 (-lmp(O) - |p(0)|) < k < 2(-Imp(0) + |p(0)|), 
D = 0, if k = 2(-Imp(0) ± |p(0)|), 

D < 0, if k < 2 (-lmp(O) - |p(0)|), or k > 2 (- lmp(0) + |p(0)|). 
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Depending on the sign of D, the nature of the semigroup ip t changes qualitatively. One 
can check that 

{hyperbolic, if 2 (-lmp(O) - |p(0)|) < k < 2(-Imp(0) + |p(0)|), 
parabolic, if k — 2 (— Imp(O) ± |p(0)|), 
elliptic, if k < 2 (-lmp(O) - \p(0)\) , or k> 2(-Imp(0) + |p(0)|). 

In the elliptic case (J7|) leads to the following equation for ipt(z) 



(8) 



\T~D + 2iAi/;t(z) + 25(1 - ip t ( z )) - k y/^D - 2%Az - 25(1 — z) + k 



yf^D - 2iAi/) t (z) - 25(1 - *l) t {z)) + k V^D + 2iAz + 25(1 — z) — k 
A necessary and sufficient condition for the curve ipt( z ) to be closed is the existence of 
a value t of t, such that ipt ( z ) — z an d V'to = ^0(2), i- e - ; both the position and velocity 
coincide with the position and velocity at the initial moment t — 0. It follows from (jSJ) 
that ?/>i(z) returns to the point z for t = m G Z. From the original equation it 

follows that the velocity of a point of the trajectory ipt{z) passing through the point z is 
^—n^- p(z I 'e lkt ) , i. e., it is ^-periodic. Thus, it is necessary and sufficient that there 
exist integers m, n G Z, such that 4=p = ^r 1 , and we arrive at the following proposition. 

Proposition 3. Let {ipt{ z )}t>o be an elliptic semigroup of automorphisms, i. e., ipt{ z ) 
solves the initial value problem, 

ipo(z) = z, 

fork<2 (-5 - VA 2 + B 2 ) , or k > 2 (-5 + VA 2 + 5 2 ). 
Then the curves ipt{ z ) ore closed if and only if 

k 



s/-AA 2 + ABk + k 2 
Figures |2] and [3] illustrate Proposition [3j 

3.2. Sufficient conditions for ellipticity. We answered the question about the values 
of k for which the semigroup {ipt{z)}t>o is elliptic, hyperbolic or parabolic completely in 
the case of automorphisms. Let us consider the general case of semigroups of holomorphic 
self-maps of D 

In the case of k = the semigroup {ipt( z )}t>o is either hyperbolic (if Zlim 2 ^i(l — 
z)p(z) < 0) or parabolic (if Zlim z _ > .i(l — z)p(z) = 0). 

In the rest of this section we assume k ^ 0. Since the semigroup {ipt( z )}t>o satisfies the 
equation 

j t U z ) = (M z ) - V 2 P(M z )) - ityt(z), 

it is necessary and sufficient for ellipticity of i/j t that the equation 

(9) (1 - z) 2 p(z) -ikz = 
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Figure 2. The closed curve t (O) for p = k = 2.5. 




Figure 3. The closed curve </>t(0) for p = i, k = 0.5. 

has a solution z inside D. 

Equation flU]) can be rewritten as 

(10) p(z) = K k (z), 

where the function 

Kk{z) = j^w 

is the rotated and rescaled Koebe function n~^ip • The i ma g e of the unit circle T under 
is the ray along the imaginary axis going from ik to oo. The inverse function is given by 
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If we rewrite fllOl) as 



K k 1 (p(z)) = z, 

then the question about zeros of the infinitesimal generator of ip t reduces to the question 
about fixed points of the map F(z) = K k 1 (p(z)) in O. 

Proposition 4. Let Kep(z) > e > for all z6D. Then ipt{z) is elliptic for all k ^ 0. 

Proof. Solutions to the equation K k 1 (p(z)) = z in the closed unit disk D must lie in 
K k \p(B)). The intersection K k \p(D)) n T may contain at most one point z = 1, and 1 
cannot be a solution of fl9]) for 7^ 0. Hence, the infinitesimal generator must have a zero 
inside D, and the semigroup {ipt(z)}t>o is elliptic. □ 

The same argument leads to a more general statement. 

Proposition 5. Let for a given k the Euclidean distance dist(p(T), Kfc(T)) > e > 0. Then 
ipt{z) is elliptic for this value of k. 

Next proposition states that, in general, large values of k correspond to elliptic semi- 
groups. 

Proposition 6. For each Herglotz function p(z) there exists k , such that for all \k\ > 
ko the semigroup {ipt{z)}t>o is elliptic. Moreover, the D enjoy- Wolff point of ipt{z) can 
approach arbitrarily close, provided that ko is large enough. 

Proof. Let B> r = {z : \z\ < r}. It is easy to see that for k large enough the map 
(11) F(z)=K k 1 (p(z)) = 



Ap(*) + 1 " 1 



is a strict contraction. Indeed, the set p(D r ) is bounded, so for large k, the set ^-p(O r ) is 

contained in an arbitrarily small neighborhood of 0, the sets ^p(B> r ) + 1 and \J ^p(D r ) + 1 

are contained in arbitrarily small neighborhoods of 1, and F(B r ) is contained in an arbi- 
trarily small disk centered at the origin. Thus, F is a strict contraction and has a unique 
fixed point Zq = F(zq), and the fixed point z lies in a small disk centered at the origin. □ 

Proposition 7. Suppose that for a given k ^ the semigroup ipt{z) defined by is 
elliptic, and let zq be its D enjoy- Wolff point. If k > 0, then Imzo — 0, and if k < 0, then 
Imzo > 0. Moreover, the equality Imz = is possible, if and only if, p(z) = Bi, Bel. 

Proof. Since the Denjoy- Wolff point of ipt{z) is the fixed point of F(z) defined by (fTTjl . then 
it must lie in the domain F(B). For p(z) = Bi, F(B) = {1 + ^ (k - Vk 2 + AB) } C R. 
Otherwise, ImF(D) > 0, and we obtain the result. □ 

Of course, it might happen that the semigroup ip t has boundary fixed points, even though 
it is elliptic. For example, let p{z) = Then the function F(z) = K k ~ 1 (p(z)) has fixed 
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points on T for all k. Indeed, the equation 

1 - e ie ike i 



which can be rewritten as 



has solutions for any k. 



l + e ie ( e ^_i)2- 



6 1 k 
tan - = -- 



2 2cos#-l 



4. Stochastic case 

In this section we consider the generalized Lowner evolution driven by a Brownian par- 
ticle on the unit circle. In other words, we study the following initial value problem. 

(12) ( **.<*, «) - «),«,«-), oo,2eD,weJ!. 

[(po(z,uj) = z, 

Here, Q is the sample space in the standard filtered probability space (Q, J 7 , (J-" 4 ),P) of 
Brownian motion, r(i, oS) = exp(ikB t (uj)), where fcel, and B t {oj) denotes the standard 
1-dimensional Brownian motion. The function p(z,t,u) is a Herglotz function for each 
fixed oj G Q. In order for <f> t (z,u) to be an Ito process adapted to the Brownian filtration, 
we require that the function p(z, t, u) is adapted to the Brownian filtration for each z G D. 

For each fixed u G Q, equation ( 1T2|) may be considered as a deterministic generalized 
Lowner equation with the Berkson-Porta data (r(-, u),p(-, •, u)). In particular, the solution 
4>t(z, u) exists, is unique for each t > and u G Q, and moreover, is a family of holomorphic 
self-maps of the unit disk. 

The equation in (|T2|) is an example of a so-called random differential equation (see, for 
instance, [27J). Since for each fixed u> G Q it may be regarded as an ordinary differential 
equation, the sample paths t h-> (j) t (z,u) have continuous first derivatives for almost all 
w (Figure H]). This is entirely different from the case of Ito and Stratonovich stochastic 
differential equations, where sample paths are nondifferentiable and tools of stochastic 
calculus are needed to solve those equations. 

4.1. An explicitly solvable example. Let p(z,t,u) = = JhB t t^)_ z -i that is, the 

stochastic counterpart of the Herglotz function in Example 1 from the previous section. It 
makes equation f|T2|) linear: 

j t Uz^) = j kBt{w) -H^), 

and a well-known formula from the theory of ordinary differential equation yields 



Jt(*,w) = e * 



z+ f e s e lkBs{uj) di 
Jo 
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Figure 4. A sample path of <f> t (0, u) for p(z, t) 
t e [0,30]. 



r(t)+z 
r(t)-z 



T(t) 



e ikB >, k — 5, 



Taking into account the fact that Ee ikBt ^ = e 2 tfc2 ; we can also write the expression for 
the mean function E(p t (z,u) 



Thus, in this example all maps <p t and E<ft t are affine transformations (compositions of a 
scaling and a translation). 

In general, solving the random differential equation ([12]) is much more complicated than 
solving its deterministic counterpart, mostly because of the fact that for almost all u the 
function t \— > r(t,u) is nowhere differentiable. 

4.2. Invariance property. If we assume again, as we did in the deterministic case, that 
the Herglotz function has the form p(z,t,u) = p(z/r(t,uj)), then it turns out that the 
process <f) t {z,u) has an important invariance property, that were crucial in development of 



(13) 




otherwise. 



SLE. 



Let s > and introduce the notation 




Then (f>t(z) is the solution to the initial-value problem 




■p[4> t (z,u)/f(t) 



) 
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where f(i) = r(s + t)/r(s) = e lk ( B s+t-B s ) j g a g a j n a Brownian motion on T (because 
B t = B s+t —B s is a standard Brownian motion). In other words, the conditional distribution 
of 4> t given <f> r , r e [0, s] is the same as the distribution of 4> t . 

4.3. Stochastic change of variables. By the complex Ito formula (Appendix [A]), the 
process = e~ lkBt satisfies the equation 

h 2 

de -ikB t = _ ike -ikB tdB e ~ ikBt dt. 

2 

Let us denote by * t ( z,oj). Applying the integration by parts formula to \E^, we 

arrive at the following initial value problem for the Ito stochastic differential equation 



(14) 



d^t = -ik%dB t + + (% - l) 2 p^te ikBt{u] \t,u)^ dt, 

®o(z) = z. 

A numerical solution to this equation for a specific choice of p(z, t, u) is shown in Figured 




Figure 5. A sample path of \l> t (0, u) for p(z, t) = i, k = 1 and t e [0, 2]. 

Analyzing the process ^j^jr instead of the original process <f) t (z,u) is in many ways 
similar to one of the approaches used in SLE theory. Recall that many properties of SLEfc, 
that is, of the solution to the problem 



9o{z) = z, 



can be proved using the fact that gt(z) — VkB t is a so-called Bessel process. 
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The image domains ty t (B,U)) differ from <p t (B,uj) only by rotation. Due to the fact that 
|\I/ t (z,u;)| = \<fi t (z,u})\, if we compare the processes <f> t (0,u) and ^ t (0,oj), we note that their 
first hit times of the circle T r with radius r < 1 coincide, i. e., 

inf{t > 0, |&(0,w)| =r} = mf{t > 0, \%(0,co)\ = r}. 

In other words, the answers to probabilistic questions about the expected time of hitting 
the circle T n the probability of exit from the disk D r = {z : \z\ < r}, etc. are the same for 
4>t(0,cj) and %(0,cj). 

4.4. ^ t as a diffusion. Associated semigroup of operators. If the Herglotz function 
has the form p(z,t,u) = p(z/r(t,u)), then the equation (fl"4"]) becomes 



(15) 




-ik%dB t + (-f * t + (% - l) 2 p(^ t )) 



dt, 



z, 



and may be regarded as an equation of a 2-dimensional time-homogeneous real diffusion 
written in complex form. This implies, in particular, that is a time-homogeneous strong 
Markov process. By construction, ^t{ z ) always stays in the unit disk. 

Let Bi,(D) denote the Banach space of bounded complex- valued Borel functions / : D — > 
C with supremum norm ||/|| = sup 2gE) 1/(^)1 . To a time-homogeneous Markov process ty t 
taking values in D one can associate a family of operators T t acting on S&(D) whose action 
is defined by the formula 

(16) (T t f)(z) = E/(* t (z)) 

(see PUD]). 

When t — 0, the operator T t is the identity map. For all t > 0, the operators T t are 
linear by linearity of expectation, and ||T t || < 1, because 

||T t /|| =sup|E/(* t (z))| < supE|/(* t («))| < sup \f(z)\ = ||/||. 
zeo zgd zeo 

Since ty t is a time- homogeneous Markov process, the operators T t form a semigroup 

mm- 



(17) T s+t = T s T t . 

The operators T t preserve the class of bounded continuous functions and the Banach 
space of bounded analytic functions if 00 (ID), as stated in the two propositions below. 

Proposition 8 (cf. [2], proof of Theorem 6.7.2]). If f : D — >■ C is a bounded continuous 
function, then T t f is bounded and continuous for every t > 0. 

Proof. It has been shown above that T t is a bounded operator. It remains to show conti- 
nuity. 

Let {z n }^ =l C D and z n — > z e D. Since f(z) and ty t (z,cu) are continuous, we have that 
f(^ t (z n ,u)) = f(^ t (z ,u)) for every u G Q. Since / is bounded we also have that 
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\f(^ t (z n ,u))\ < ll/H G L 1 (r2,(iP). Hence, by the dominated convergence theorem, 
hm(T t /)(z n ) = hm Ef(%{z n ,u)) = Ef(%(z,u)) = (T t f)(z ) 

□ 

Proposition 9. // / G H°°(B) then T t f G H°°(z) for every t > 0. 

Proof. Let / : D — > C be a bounded analytic function. T t f is continuous, and thus, 
by Morera's theorem, it suffices to show that f. Ef($f t (z), oS) dz = for every triangle 
A C D. But this follows directly from Fubini's theorem. Indeed, let 7 : [0, 1] — > B be a 
parameterization of the triangle A. Then, since 

E / |/(**(7W,w))Y(0l*< \\f\\-l(l)< 00, 
Jo 

we conclude that 



/ Ef(%{z,u))dz = E f f{%{z,u))dz = 0. 

J A J A 



□ 



4.5. Infinitesimal generator of ty t - The operator A defined by the expression 

Af( , y (T t f)( Z ) ~ f( Z ) 

Af(z)=\un , zeB, 

acting on a suitable class of functions / : C — > C is called the infinitesimal generator of 
the Ito diffusion ty t . The set of functions for which Af is defined is denoted by T>a, i- e. 
/ G T>a if the limit above exists for all z6D. 

Following the lines of the standard proof for the case of real diffusions and applying 
the complex Ito formula (see Appendix [A} one shows that Cq(C) C T>a, and that the 
infinitesimal generator in this case is given by the formula 

A = (4« + (,-!)•«,)) 1 

° L 1,2=2 ° , 7-21 _|2 a J 



+ [--z+(z- iyp(z) - - -k^— + k*\ z 



2 J dz 2 dz 2 dzdz 

Now, let / be analytic in D. Then, by Ito's formula for analytic functions, 



E/(* t (z)) - /(*) 

t r / p \ « 1,2 o2 x 



= E 

and thus, 

(18) a=(-^ + (,- l»)Ji;-i^V 



tit, 
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Let p(z) = Yl^=o a n zTl be the Taylor expansion of p(z) in D. Set 

£„ = -z n+1 |-, ra = 0,±l,±2,... 

oz 

Observe that C n give the Virasoro generators in a representation of the Virasoro algebra 
over the space of complex Laurent polynomials. Interesting connections between highest- 
weight representations of the Virasoro algebra and properties of SLE were studied, e. g., 
in [31 12] 

Using the formula z 2 = + £ , we can rewrite f fT8l) as 
(19) A = - 22 ( a n+i ~ 2a n + a n _i)£ n - y^o, a -i = a -z = °- 

n=-l 

Note that operators C n with n < — 1 never appear in this expression for A. 

We can now write the generators for particular functions p(z) considered in the previous 
section. 



For p(z) 
For p(z) 
For p(z) 



l 

1-2 > 



1+2 
1-2' 



A = + £ - y 



2 



A; 2 

A = — iC-i + 2i£ — i£i — 2~^o- 



Since the first term in (Tl9|) depends linearly on and the second term does not depend 
on it, it is easy to obtain the generator for p(z) = A + Bi: 

k 2 

A = -(A + Bi) £_! + 2iB£ + {A- Bi) C x - —C 2 . 

4.6. Analytic properties of the family T t f. Kolmogorov's backward equation. 

Let / be a bounded analytic function, / e H°°(3). From Dynkin's formula 

E[f(%(z))} = f(z) + E \ f{Af){*.{z))d8 

Jo 

it follows that the function u(t,z) = (T t f)(z) is different iable with respect to t, and 

£u(t,z) = E[{Af){* t (z))). 

The family of functions (T t f)(z) is a uniformly bounded family of holomorphic functions 
in D. Hence, it is normal by Montel's theorem. Since this family is pointwise continuous, 
it is also continuous in the topology of local uniform convergence by Vitali's theorem. 

The function u(t,z) = (T t f)(z) satisfies the following initial value problem for the Kol- 
mogorov backward equation. 
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fz+(z-l) 2 p(z)) 



du _ 1l2 2d 2 u 

dz 2^ ^ ai 2 "' 



z e D. 



4.7. Relation to the deterministic case. Let us denote by u(t, z) the mean function 
of ^t{ z )-, which is the result of acting by T t on the identity function: u(t, z) = (T t id)(z) = 
E^f t (z). It follows from ([TTJ) that 



This reveals a deep similarity between the deterministic and stochastic cases. While the 
deterministic family i[>t(z) of holomorphic self- maps of the disk forms a bona fide semigroup, 
the stochastic family \l/ t forms a semigroup only in average. 

We can also look at this connection from a different point of view: to the semigroup 
ipt{z), as well as to the diffusion ty t , one can associate a family of operators T t , which in 
the deterministic case are defined trivially by (T t f)(z) = f(tp t (z)). 

4.8. Zero points of the stochastic vector field. Analogously to the deterministic case, 
we can regard ( I14p as an autonomous Lowner equation 



Here, W t (u) denotes a generalized stochastic process known as white noise (see [18] for 
details). 

By the argument we used in subsection 13.21 it follows that for k ^ there exists a unique 
point Zo G D, such that — \z§ + (zq — l) 2 p{zo). 

By the definition of the infinitesimal generator of a diffusion, zq is the unique point in 
the unit disk having the property 



4.9. Hierarchy of differential equations for moments. We use the term the n-th 
moment of the process ^t{z) for the function fi m (z) = T t z m = E\l/™(z). The first moment 
of is thus, the mean function Ety t {z)- 

By applying the Ito formula to the function f(z) = z m and taking expectations we obtain 



For a function p(z) = Y2^=o a nZ n this leads to an infinite-dimensional system of differ- 
ential equations for the functions /ii, // 2 , ■ ■ ■ : 



u(t + s,z) = Eu(Q, Vt+s{z)) = Eu(t, $> s {z)) = Eu(s, %(z)). 




lim E 



t 



jz + (z-l) 2 p(z). 
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(20) — /2 m = a mij, m _i + (a,! -2a -m y J m(i m + ^(a n _i - 2a n + a n+ i)mfi. 

^ ' n=l 

In the simplest case when p(z) — this reduces to 



m+n- 



d ( k 2 , , 

—fJLm = mum-i - m \ m Y + 1 ) / im ' 



and for m — 1 can be easily solved explicitly: 

1 _ e -(l+fc 2 /2)t 



As a by-product, we can find the expression for the covariance of the processes f Q e s e lkBs ds 
and e~ lkBs simply by comparing the formula above with the formula f|T3|) for E<p t (z) : 

Ee ~ikB t = e -4\ 



e s e~ ikB 'ds 



e~% k 2 = 2 



!•->/ 



l-k 2 /2 ' 



otherwise, 



t \ i _ -(l+fc 2 /2)i 

e s e-^cfe • e" ifcBi ) = — . 

; l + ^ 2 /2 



Cov 



e s e lkB °ds,e 



ds ■ e 



1 l„-3t 



e" 3 *(3t + 1), 



A; 2 = 2, 



l_p-(l+fe /2)i -tfe / /2_ -t _1«,2 . 

l+ifeV2 1^72- e 2 ' 0th ™ Se - 



Even though it is impossible to solve ( 12"U|) explicitly for other Herglotz functions p(z), the 
system still allows in many cases to express higher-order moments in terms of lower-order 
moments. If, for instance, p(z) = j^, then the system (1201) is written 



and hence, 



d 
dt 



k 2 m 2 



k 2 m 1 d 

2 mat 
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4.10. Growth estimates. First, we write the equation ([12]) in polar coordinates using 
4 t |, 9 = arg4>t, U = Rep, V = Imp : 

- kB t )(l + r 2 ) - 2r] U(r, 9) + sm(9 - kB t )(l - r 2 )V(r, 9), 
(i - r) sin(# - kB t )U(r, 9) + [(1 + r) cos(fl - fcS t ) - 2] K(r, 0). 

Let the Herglotz function again be given by p(z) = i±f Rep(0) +ilmp(0). Then (]2T|) 
becomes 




(22) 



| = (1 _ r 2 )\p(0)\ cos(# - kB t - axgp(0)), 



§ = -2Imp(0) - (i + r) |p(O)|siii(0-A;S t -argp(O)). 
Let r = 100(^)1 — M- Then the solution to the first equation in (|22|) is given by 

r(£) = tanh ^|p(0)| f*cos(9 — kB s — argp(O)) ds + arctanhr j , r < 1, 



(23) 



r(i) = 1, r = 1. 



In particular, this implies that (j)t(z) maps T onto T, and the function p(z) = Rep(O) + 
ilmp(O) generates families of automorphisms of D in the stochastic case, too. 
From ( 123]) it is easy to obtain the following estimate 

(24) tanh(— |p(0)| t + arctanhr ) < r(t) < tanh(|p(0)| t + arctanhr ). 

Note that the estimate from below quickly becomes trivial as t increases (tanh takes 
negative values, but r(t) always stays non- negative). 

In a similar way we can obtain estimates of growth of r(t) for other Herglotz functions. 
For p(z) = 

(25) r e~ l - (1 - e _t ) < r(t) < r e~ l + (1 - e _t ). 
For p[z) = 1, 

(26) < r(t) < r ° (1 7 f) + '. 
^ ' l+t(l + r„) " w " l+t(l-r„) 

It can be checked that the estimates from above in ( 124]) (in the case p(z) 
and ( 126]) are sharp, and the equality is attained when k = 0. 



l±z\ 



4.11. Boundary diffusion. Suppose again that the family <fit{z) is a family of disk au- 
tomorphisms, p(z) = Rep(O) + lmp(0)i. Since the unit circle T is always mapped by 
4>t{z) bijectively onto itself, we can study the stochastic family of diffeomorphisms of the 
unit circle described by the equation 

d9 

— = -2Imp(O)-2|p(O)|sin(0-A;B t -axgp(O)), 

which is the restriction of the second equation in (122]) to T. 

The process Qt = arg ~ argp(O) = 9 — k B t — argp(O) is a diffusion on T = M/27tZ: 

(27) dB t = -2 (Imp(O) + \p(0) | sin Q t ) dt - kdB t , 
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with the infinitesimal generator given by 

A = -2 (Imp(O) + |p(0)| sin#) — + — dB t . 

The functions annihilating the infinitesimal generator, that is, the solutions to the equation 
Af = are given by the formula 

f(9)=c 1 + c 2 J exp — ds 

The function p(z) = corresponds to the noisy North-South flow 

(28) dB t = -2 sin Q t dt - kdB t , 
studied in E] (see also [21] 125]). 

Appendix A. Complex Ito formula and integration by parts 

We recall two basic results of stochastic calculus in the complex plane. 
Let Z t = {X}, Yf 1 , . . . , X™, Y™) be a 2n-dimensional real Ito process, which may also be 
regarded as an n-dimensional complex Ito process, by identifying Z\ = Xf + iY t k . Let 

f(z x , ...,z n ) = (u^zx, ...,z n ), Vtfa, . . . [Zi, . . . , Z n ),V m {Zi, . . . , z n )) 

be a C 2 function from M 2n to M 2m . Then it follows from the real multidimensional Ito 
formula that W t = f(Z t ) is a 2m-dimensional Ito process, and its p-th (complex) component 
is given by 

(29) dW p = d(Re W p + i ■ Im W p ) = jr^ dx j + Yl IT dY i 

j=i j j=i "Vi 

+ lT,lf^ dx rdX k + Y 1 f^dX r dY k + lj2p^ dY rdY k 

2 ^ dxjdxk dxjdy k 2 dyjdy k 

3->k 3}to 

1 dx , x y x d Y k + ly J^dYj ■ dY k ] , 

2^dx,dx k 3 ^d Xj dy k 3 2^dy j dy k 3 k ) 

0,k J 3,k J j,k / 

where the dot product of stochastic differentials is a linear operation computed according 
to the following formal rules dBi ■ dBj = dijdt, dBi ■ dt = dt ■ dBi = 0. 

If we use the notations |£ = § ( |^ - i-jjt) and = ± (-^- + i^f-) , then (1291) sim- 

dz k 2 \dx k dy k J dz k 2 \dx k dy k J ' V " 

plifies to the well-known complex Ito formula (see, e. g., |28j): 
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1 n r) 2 f 1 n f) 2 f - r) 2 f 

j,k J j,k J j,k J 

Let Z t , W t be two Ito processes in the complex plane. Take f(z, w) — z ■ w. Since all 
terms containing J| and ^ vanish, we obtain the integration by parts formula for complex 
Ito processes: 

d(Z t W t ) = W t dZ t + Z t dW t + dZ t ■ dW t . 

In general, if the process Z t takes values in a domain DcC" where the function / is 
holomorphic, the formula (|30|) reduces to 

3=1 J j,k J 
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